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ABSTRACT 



Context. The microquasar Cygnus X-3 was detected at high energies by the gamma-ray space telescopes AGILE and Fermi. The 
gamma-ray emission is transient, modulated with the orbital period and seems related to major radio flares, i.e. to the relativistic jet. 
The GeV gamma-ray flux can be substantially attenuated by internal absorption with the ambient X-rays. 

Aims. In this study, we examine quantitatively the effect of pair production in Cygnus X-3 and put constraints on the location of the 
gamma-ray source. 

Methods. Cygnus X-3 exhibits complex temporal and spectral patterns in X-rays. During gamma-ray flares, the X-ray emission can 
be approximated by a bright disk black body component and a non-thermal tail extending in hard X-rays, possibly related to a corona 
above the disk. We calculate numerically the exact optical depth for gamma rays above a standard accretion disk. Emission and 
absorption in the corona are also investigated. 

Results. GeV gamma rays are significantly absorbed by soft X-rays emitted from the inner parts of the accretion disk. The absorption 
pattern is complex and anisotropic. Isotropization of X-rays due to Thomson scattering in the companion star wind tends to increase 
the gamma-ray opacity. Gamma rays from the corona suffer from strong absorption by photons from the disk and cannot explain the 
observed high-energy emission, unless the corona is unrealistically extended. 

Conclusions. The lack of absorption feature in the GeV emission indicates that high-energy gamma rays should be located at a 
minimum distance ~ 10^-10'" cm from the compact object. The gamma-ray emission is unlikely to have a coronal origin. 



Key words, radiation mechanisms: non-thermal - 
disks 
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1. Introduction 

Gamma-ray emission from Cygnus X-3 was detected by the 
gamm a-ray space telescopes AGILE and Fermi ( Tava ni et al.l 
120091; lAbdo et al. 200 9|). This syste m is formed of a Wolf- 
Rayet star (van Kerkwiik et alj|1996h and an unidentified com- 
pact object (a black-hole or a n eutron star) in a very tight 
4.8 h orbit, situated at 7-9 kpc (iDickevI 119831: iPredehl et al] 
l2000HLing et al. 20091 Cy gnus X-3 is also known as a micro- 
quas ar as there is clear evidence of resolved re lativistic radio 
jets jMarti et al.ll2000l; iMioduszewski et alj|2001 l). The gamma- 
ray emissio n is transient, mo dulated with the orbital period of 
the system (Abdo et al. 2009) and seems to be related to giant 
radio flares and episodes of major ejections of a relativistic jet. 
The spectrum measured by Fermi is a soft power-law of index 
~ 2.7 with no sign of high energy cut-off, for a total luminos- 
ity above 100 MeV of about 3 x 10^* erg s^' (at 7 kpc). The 
gamma-ray orbital modulation was interpreted as the result of 
Compton upscattering of thermal photons fro m the Wolf-Rayet 
star by relativistic pairs accelerated in the jet (lAbdo et al .1120091: 
iDubus et al.l l2010'). 

The propagation of energetic gamma rays to the observer 
can be greatly affected by pair production j + j e^ + 
e^ with the ambient radiation field in the system, regard- 



less the underlying production mechanisms for these gamma 
rays. This issue motivated many different studies in the past, 
at a time when Cygnus X-3 was thought to emit GeV, TeV 
and even PeV gamma-ray photons (see for instance IVestrandl 
1983L iGouldl 1983; Ancarao 1984; Cawlev & Weekes 1981 



Protheroe & Stanev 1987t lCarramifianalll992l: IWu et al.lll993l) . 
Hermsen et al. ( 1987) invoked the effect of gamma-ray absorp- 
tion as a possible explanation for the non-detection of Cygnus 
X-3 at GeV energies by COS-B. This question is quantitatively 
revisited here in the context of the new gamma-ray observations 
of Cygnus X-3. The goal of this article is to formulate constraints 
on the location of the gamma-ray source in the system. 

High-energy gamma rays in the energy band ei = 100 MeV- 
10 GeV annihilate with soft X-rays photons of energy emin ~ 
m^c'^/ei = 25 eV - 2.5 keV. Hence, the effect of the gamma- 
ray opacity at high energy depends on the ambient X-ray density 
in the system. Known Wolf-Rayet stars have effectiv e temper- 
atures around ~ 10^ K (see e.g. the review by ICrowtherl 
|2007|) and provide a high density of ~ 20 eV < emin thermal 
radiation (with a bolometric luminosity > lO-'^ erg s"'). 
Absorption with stellar radiation is then irrelevant in this con- 
text. Note that the star would have a major role in the op- 
tical depth for very-high energy (> 100 GeV) gamma rays 
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(lFordlll984t iProtheroe & Stanevlll987t iMoskalenko etanil993t 
lBednarekll2010h . Some Wolf-Rayet stars present an excess in 



soft X-rays with luminositi es up to ~ 10 erg s" , particularly 
in nitrogen-rich stars (e.g. Skinner et al.ll2010h . This source of 
photons would be negligible in Cygnus X-3 where the X-ray 
luminosity is about 10^^ erg s~\ and will not be considered 
here. If the accreting compact object is a neutron star, thermal 
radiation is emitted by the surface of the star With a typical 
surfac e temperature T^s ~ 10* K (see e.g. lYakovlev & Pethickl 
l2004l) . the neutron star provides ~ 0.25 keV soft X-rays able 
to absorb 1 GeV photons. The luminosity of the neutron star 
Lns ~ 10^^ erg s"' is too small and will be neglected as well in 
the following. The dominant source of soft X-rays in Cygnus X- 
3 is provided by the accretion disk formed around the compact 
object. The inner parts of the disk could be very hot with temper- 
ature Tdisk ~ 10^-10^ K and emit copious ~ 1 keV soft X-rays in 
the system. We provide in Sect.|2]and in AppendixlAla complete 
study of the absorption of GeV photons in the radiation field of 
an accretion disk. 

Long-term observations in X-rays revealed tha t Cygnus X-3 
presen ts complex spectral and temporal features. ISzostek et alj 
(l2008h classified the X-ray spectral states into five distinct 
groups, from the 'hard state' to the 'soft state'. In the hard state, 
the X-ray flux is dominated by hard X-rays (> 10 keV) and 
the spectrum is well-fitted by comptonized emission on ther- 
mal and non-thermal electrons with a cut-off at « 20 keV and 
a Compton reflection component. The soft state is dominated by 
the bright thermal emission from the accretion disk in soft X- 
rays (~ 1 keV), and contains also a non-thermal tail in hard X- 
rays extending at least up to 100 keV. Gamma-ray flares were ob- 
served only when Cygnus X-3 was in the soft state (during major 
radio flares). Therefore, we focus only on the soft state and ex- 
amine whether the non-thermal X-ray component could, in addi- 
tion to the accretion disk, contribute significantly to increase the 
gamma-ray opacity in Cygnus X-3 (Sect.O. The nature and the 
origin of this non-thermal emission is unknown. This is possibly 
related to electrons in the hot corona above the disk or in the rel- 
ativistic jet. The extrapolation of the hard X-ray tail in the soft 
state up to the GeV energy domain is tempting. In Sect. |3] we 
examine also the possibility that the gamma-ray emission comes 
from a hot corona above the accretion disk using the one-zone ra- 
diative model BELM described in Belmont et al. (2008). Gamma- 
ray emission_Jrom a magnetized corona was also investigated 
by iRomero et al.l (f2010.) with an application to the microquasar 
Cygnus X- 1 . 

The global and simplified picture of Cygnus X-3 adopted in 
this paper is summarized and depicted in Fig. [1] The main re- 
sults and conclusions of this study are exposed in the last section 
(Sect.©. 



2. Absorption by the accretion disk 

In this section, we propose a full study of gamma-ray absorption 
by the photons from an accretion disk formed around the com- 
pact object (Fig. [I])- With simple assumptions on the physics and 
geometry of the disk ( 32.11 ). we derive the main properties of 
the optical depth for gamma rays ( 32.211231 ) with an application 
to Cygnus X-3. Soft X-rays emitted by the disk can be signifi- 
cantly scattered and absorbed by the dense Wolf-Rayet star wind 
in Cygnus X-3. We discuss in the last section ( 32.41 ) the impact 
of the stellar wind on the gamma-ray opacity. 



observer 




counter 
jet 



Fig. 1. Diagram of the microquasar Cygnus X-3 with the dif- 
ferent sources of photons considered in this study. The ob- 
server sees the disk at an angle t//. In cylindrical coordinates, 
the gamma-ray source is located at {r,(f>s,z), possibly in the rel- 
ativistic jet. For a close-up view of the disk see the appendix, 
Fig.lAH 



2. 1 . Accretion disk model 

Pair production of high-energy gamma rays in the radiation 
field of an accretion disk has been investigated by several 
authors in the c ontext of X-ray binaries (Carraminana 119921: 
Bednarek'1993'), but also in AGN (see e.g. Becker & KafatosI 
1995; Zhang & Cheng 1997; Sitarek & Bednarek 2008). In most 
of these previous studies, the accretion disk is optically thick, 
geometrically thi n, flat and follows the stand ard quasi-Keplerian 
disk-a model of IShakura & SunvaevI (Il973h . According to this 
theory, each surface element of the disk is at thermal equilib- 
rium. In the non-relativistic case and if the boundary conditions 
of the disk are neglected, the radial profile of the effective tem- 
perature follows 



3/4 



(1) 



where R = R/rg. R is the distance from the center of the disk and 
Tg = GMco/c^ is the gravitational radius of the compact object. 
Tg is formally the temperature of the disk at R = r„ which de- 
pends on the accretion rate Mace onto the central compact object 
of mass Mco', T„ is defined as 



/ 3MaccC* 



1/4 



(2) 



where G is the gravitational constant and osb is the Stefan- 
Boltzmann constant. Eq. ([1]) can be rewritten as 



T{R) = 7.7 X 10" Mg^'^M';^^/^^'* K, 



(3) 



where Mg = M^^JIO'^ Mq yr"\ M20 = Mco/20 Mq and 
^7 - R/IQ^ cm. The inner radius of the accretion disk 
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T at 1 GeV {on axis) 




Fig. 2. Optical depth Tyy map for a gamma-ray photon of energy 
ei = 1 GeV, integrated along the line of sight up to the observer 
situated at infinity. The compact object is a M^o - 20 M© black 
hole accreting at Mace = 10"** Mq yr^K The inner radius of the 
disk is Rin = 6 rg 1.8 x 10^ cm and the outer radius is /?out - 
8 X 10^ cm. The altitudes z = Rin and z = Rom are indicated 
by dotted lines. Absorption is calculated with the temperature 
profile given in Eq. (|4]i. Primary gamma rays are injected along 
the axis of the disk at an altitude z(r = 0). The observer sees the 
disk with an inclination tfr. Bright regions indicate low opacity 
or 'UNABSORBED' regions, and dark regions high opacity or 
'ABSORBED' regions. Level lines T^y = 10 ^ 10 ', 1, 10, 10^ 
are overplotted. 



is not well constrained but is usually set at the last stable or- 
bit (i.e. Rin = 6 rg for a non-rotating black hole) in the soft- 
state where the accretion disk component is very bright (e.g. 
lEsin et aT ] ll997h . For an accreting neutron star, we will assume 
that the inner part of the disk reaches the neutron star surface 
Ri„ = 7?ns = 10 km. The emission produced by the dissipation 
of the kinetic energy of the accreted material in the neutron star 
surface is not considered here. We will do the approximation that 
the observed X-ray luminosity comes from the disk only. 

The temperature profile can be substantially modified close 
to the compact object by the inner boundary condition of 
the disk, neglected in Eq. ([T]l- Using a torque-free inner 
boundar y condition and a p seudo-Newtonian gravitational po- 
tential fPaczvriski & Wiital 119801) . Eq. ([TJ is changed into 
(.Gierliriski et al.. l999.) 



t{r)^t. 



R-2/3 



r[R-2) 



1 - 



2^ ^3/2 



1/4 



(4) 



For ^ < 10, the temperature becomes significantly smaller than 
in Eq. ([TJ (about a factor 2 at ^ = 10). Far from the last sta- 
ble orbit ^ » 6, this solution matches the profile T oc R^^l'^. 
Since the effect of the gamma-ray absorption depends critically 
on the temperature in the inner parts of the disk, we compute 



here the optical depth for high-energy gamma rays in the frame- 
work of a type disk with the modified profile temperature given 
by Eq. (|4]i. General relativistic effects such as Doppler beaming 
and the bending of light trajectories are neglected. 

In close binaries, the accretion disk does not extend beyond 
a fraction of the Roche lobe radius due to tidal torques exerted 
by the companion star, (Paczvriski' 19771: iPapaloizou & Pringld 
,1977). We set the outer radius of the disk i^out at the Roche lobe 
radius. The energy of the photons emitted in the outer part of 
the disk falls in the infra-red band, hence is too small to absorb 
100 MeV - 10 GeV gamma rays, but could annihilate with TeV 
photons. 

The total luminosity of the disk Ljisk depends on the accre- 
tion rate of matter Mace onto the compact object, such that (for a 
disk truncated at < R extending at infinity and following the 
temperature profile in Eq. [1] Longair 1994) 



idisk = — — ~ 8.4 X 10" MsM2oRj' erg s 



2Ru 



(5) 



We assume for simplicity that the total bolometric luminosity ob- 
served in X-rays (corrected from absorption by the Wolf-Rayet 
star wind and the interstellar medium) originates from the disk 
in Cygnus X-3, so that Lx ~ idisk- Pairs produced by the ab- 
sorption of gamma rays could interact with photons from the 
disk or from the star and develop an electromagnetic cascade, 
reduc ing the opacity in th e system (s ee Akharonian & VardanianI 
19851 ICarramifiana([T99l IWu et al.l[T993; ■Bednarek 2010) . The 
cascade might be quenched if the ambient magnetic field is 
large enough in the system (see Sect. 13.2b . Synchrotron emis- 
sion from secondary pairs may be i mportant or even dominant 
in hard X-rays and soft gamma rays (Bosch-Ramon et al.ll20()8t 
ICerutti et al.ll2010tlRomero et al.ll2010 l). The cascade above the 
accretion disk is not considered in the following. 



2.2. Anisotropic effects 

The radiation field from the accretion disk seen by a gamma ray 
is highly anisotropic. The cross section for pair production de- 
pends on the scattering angle between the gamma ray and the 
target photon. The gamma ray is injected above the accretion 
disk at a location given by (r,<p^,z) in cylindrical coordinates 
(see Fig. lA.lb . We compute the total gamma-ray optical depth 
exactly, i.e. by including all the geometrical effects. AppendixlAl 
provides the full details for the calculation of the gamma-ray op- 
tical depth Tyy. Fig.|2]is the result of the numerical integrations 
in Eq. ( lA.lb . and gives the total optical depth above the accre- 
tion disk as a function of the viewing angle ^ and the altitude z 
where gamma rays are injected. The temperature profile of the 
disk follows Eq. (HJi. In this calculation, the gamma-ray source 
is on the axis of the disk (r = 0, see Figs. [1] lA.lb and emits 
1 GeV photons. For illustrative purpose, we consider the case 
of a 20 Mq black hole accreting at a rate Mace - 10"^ Mq yr"' 
(corresponding to Ldisk ~ 10^^ erg s '). 

The altitude at which the medium becomes transparent to 
gamma rays is smallest at low inclinations < 30°). The direc- 
tion of the gamma rays and the photons from the disk are about 
the same for these inclinations. For higher inclination {ijj > 30°), 
the GeV photon crosses a thicker radiation field and the colli- 
sion with the photons from the disk becomes closer to head-on, 
resulting in a higher opacity for a given altitude z. In this case, 
gamma rays should be injected at z > 10**-10^ cm for = 90°. 

The gamma-ray optical depth depends also highly on the en- 
ergy of the gamma-ray photon. The map shown in Fig.|2]changes 
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Black hole (^=30°) Neutron star {i/=70°) 
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Fig. 4. Optical depth map for 1 GeV photons in Cygnus X-3, as a function of the location of the gamma-ray source above the 
accretion disk of coordinates (r,0s,z). The level lines corresponding to Tyy - 10"^, 10 1, 10, 10^ are also overplotted. The 
accreting compact object is a 20 M© black hole (i.e. iff = 30°, left panel) or a neutron star (t// = 70°, right panel), located at the center 
of the disk (origin). The dashed line labelled 'ISO' shows the contour Tyy = 1 in the case where the radiation from the accretion 
disk is fully isotropized. The maps shown here are calculated for <p^ = 0° (see Fig.lTJ. The inner radius of the disk is Ri^ = 6 rg and 
the outer radius chosen at the Roche lobe radius for the compact object. Dotted lines indicate the inner and outer radius of the disk. 




e, (eV) 

Fig. 3. Transmitted flux exp(-Tyy) as a function of the gamma- 
ray photon energy ei , for the following viewing angle if/ - Q° 
(right), 30°, 60° and 90° (left). The disk has the same phys- 
ical properties as in Fig. |2] Gamma-rays are injected at z = 
100 Rin for illustrative purpose only. The Fermi-hPC[ energy 
band (100 MeV-10 GeV) is delimited by the dotted lines. 



significantly with energy, although the same anisotropic pattern 
remains. Fig.|3]shows this dependence for various inclinations. It 
can be noticed here that in addition to the massive star, the accre- 
tion disk (the outer parts emitting in optical-infrared) can also be 



a strong source of gamma-ray absorption for Te V gamma-rays 
(if an y, see the upper-limit and the discussion in lAleksic et al.l 
I2010h . 

2.3. Map of the gamma-ray optical depth in Cygnus X-3 

The inclination of the system with respect to the observer is un- 
known in Cygnus X-3. W e consider the two extreme solutions 
chosen by ISzostek & Zdziarski (.2008,) . based on their estimate 
of the minimum and maximum mass loss rate of the Wolf-Rayet 
star In the first solution the inclination of the system is / = 70°, 
corresponding to the case where the compact object is a neu- 
ti-on star (M^o -1.4 Mq) orbiting a 5 Mq Wolf-Rayet star At 
the other end, the inclination is / - 30° and corresponds to a 
20 Mq black hole and a 50 Mq Wolf-Rayet star companion. 
These two cases represent also two extreme solutions in terms 
of gamma-ray absorption as the optical depth is higher for high 
inclinations (Fig.|2]i. In addition, these two cases correspond re- 
spectively to the maximum and minimum value for the inner 
radius of the accretion disk, hence for a hotter and cooler accre- 
tion disk. Assuming that the disk is in the orbital plane (so that 
/ = iff), we investigate the spatial distribution of the optical depth 
for GeV photons above the accretion disk. 

Fig. m presents the result of the calculation, for the partic- 
ular case 0s - 0° (Figs. [H lA.ll ). We take the bolometric X- 
ray luminosity of Cygnus X-3 as the total luminosity of the 
accretion disk L disk ~ Lx, with Lx ~ 10^** erg s"' (see e.g. 
IVilhu et alJl2009l) which is about the Eddington luminosity Lead 
for the neutron star (accretion rate of Mace = 10"*^ Mq yr"'). 
In these maps, the line Tyy = 1 delimits the contour of the 
'gamma-ray photosphere' above the disk, i.e the surface be- 
yond which gamma rays can escape from the system and reach 
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Fig. 5. Azimuth {left panel) and energy dependences {right panel) of the gamma-ray photosphere contour (r^y = 1) above the 
accretion disk, for a 20 Mq black hole = 30°). In the left panel, 0s = 0°, 60°, 120°, and 180° with ei = 1 GeV. In the right panel, 
the energy of the gamma-ray is ei = 1, 10 and 100 GeV for - 0° ■ Pair production with stellar radiation is not included in this 
plot. The dotted line shows the inner radius of the disk. 



the obser ver. Note that similar maps were found in the context 
of AGN (ISitarek & Bednarekll2008l l2010l) . Absorption is more 
marked in the neutron case since the inner radius of the disk 
is Mbh/Mns ~ 14 times smaller than in the black hole solu- 
tion. Fig. |5](/e/f panel) shows the variations of the gamma-ray 
photosphere location with the azimuth angle (p^. Gamma-ray ab- 
sorption is stronger for (p^ - 180° than for 0s = 0° because the 
gamma ray photon crosses a thicker radiation field from the disk 
before it reaches the observer The gamma-ray photosphere is 
located at a distance between ~ 10"*-10" cm from the compact 
object for the neutron star and ~ 5 x 10^-10'' cm for the black 
hole at 1 GeV. Because GeV gamma rays seem not absorbed in 
Cygnus X-3, the gamma-ray source should lie beyond this criti- 
cal distance. The gamma-ray photosphere contour depends also 
on the energy of the gamma-ray photon (Fig. |5] right panel). 

The transparency constraint depends sensitively on the value 
of the inner radius of the thin accretion disk, since GeV gamma 
rays are mostly absorbed by photons emitted in the innermost, 
hottest region. For thin disks with higher gamma-ray ab- 
sorption is lower and is negligible at 1 GeV if 



Rin > 1.4 X 10^ mI'^mK^^ cm. 



(6) 



i.e. for disks with such inner radii, the mean energy of produced 
photons equals the minimal energy required for pair production 
2.7 kT = 0.25 keV (see Eq. |3] for » 6 and Eq. |AJli- 
The (unspecified) innermost accretion flow might still contribute 
to gamma-ray attenuation. However, the high X-ray luminos- 
ity would be difficult to account for if the thin disk is truncated 
above k, 50 r^. The accretion rate would have to be as high as 
about half the Eddington accretion rate limit for a 20 Mq black 
hole (MEdd = ^Edd/c^)- 



2.4. X-ray scattering and absorption in tiie stellar wind 

Thomson scattering and absorption of soft X-rays (< 10 keV) in 
the dense Wolf-Rayet star wind is a major issue in Cygnus X-3. 
This effect is partly responsible for our poor knowledge of the 
physical properties of the disk such as the temperature and the 
instrinsi c luminosity, in spit e of man y years of observations in 
X-rays dSzostek & Zdziarski 2008; Hialmarsdotter et al.ll2008h . 
The stellar wind could significantly deplete and redistribute the 
soft X-ray density seen by the GeV gamma ray, relaxing the con- 
straints formulated in the previous section (Sect. 12.3) on the min- 
imum distance to the accreting object. 

We calculate the conservative photoelectric optical depth to 
estimate the absorption of soft X-rays in the stellar wind. For a 
neutral, isotropic and radial stellar wind assumed composed of 
pure helium, the density of helium atoms at a distance A from 
the center of the star is 



"He (A) = 



47rmHyUiA2vw 



(7) 



Mw is the mass loss rate of the star, me is the hydrogen atom 
mass, Vw the stellar wind velocity and yui is the mean molecular 
weight equals to 4 for an helium wind. With a typical mass loss 
rate = 10"^ Mq yr"' and a wind velocity v„ = 1500 km s"' 
in Wolf-Rayet stars, nee ~ 5 x 10" cm"-' at the compact object 
location. The column density of target atoms A^He crossed by X- 
rays from the compact object up to a distance p is 



Nne (P) 



Jo 



nHe{A)dp', 



(8) 



where p' is the path length integration variable between the 
source of X-rays (p' = 0) and the distance p. The distance A is 
related to p', the orbital separation d and the angle 6 between the 
X-ray photon direction and the orbital plane (see Fig. lA.U via 
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the relation - p'^ + cP- 2p'd sin 0. At 1 ke V, the photoelectric 
absorption cross section is about crpE(l keV) a; 2 x 10"^^ cm^ 
(Morrison & McCammon 1983) (assuming that the abundance 
of metals per unit mass in the wind is similar to cosmic abun- 
dances). The photoelectric optical depth in the wind is T„(p) - 
A^h(p)o-pe(1 keV) > 1 if p > 10^° cm (for the black hole and the 
neutron star case). This critical distance is much bigger than the 
gamma-ray photosphere radius in the black hole case but is com- 
parable for the neutron star solution. This is probably an overes- 
timate in Cygnus X-3 as the wind should be fully ionized around 
the X-ray source over distances comparable to the orbital separa- 
tion ~ 10" cm (Szostek & Zdziarski"2008), hence much larger 
than the extension of the gamma-ray photosphere. We conclude 
that absorption of soft X-rays in the wind should not affect the 
gamma-ray opacity above the accretion disk. 

Thomson scattering between X-rays and unbound electrons 
in the wind could result in a redistribution and isotropization 
of the directions of the photons. For a fully ionized wind, the 
Thomson optical depth in Cygnus X-3 gives tj ~ 2A^HeO"T ~ 
0.5-1.8, for tfr = 30° and 70° respectively (where the factor 2 
accounts for 2 electrons per He atoms, crj is the Thomson cross 
section, and A^He is the helium column density integrated up to 
infinity in Eq. [8]l. Gamma rays would see a different angular 
distribution (more isotropic) of the incoming X-ray distribution. 
This might increase the gamma-ray optical depth (as collisions 
between photons would be on average closer to head-on) and 
reduce the anisotropic pattern of the gamma-ray photosphere 
in Fig. 3] in particular in the neutron star case. To address this 
question more quantitatively, we performed the same calcula- 
tion of the gamma-ray optical depth as in the previous section 
assuming that the angular distribution of the photons from each 
unit surface of the disk is isotropic (the integrand in Eq. lA. II 
is averaged over a uniform distribution of pitch angle do). This 
extreme situation of full isotropization yields an upper-limit to 
the effect of X-ray scattering in the stellar wind on the gamma- 
ray optical depth. Calculations show that the gamma-ray photo- 
sphere is nearly spherical of radius ^ 10'" cm (Fig. |4] dashed 
lines). The photosphere is about a hundred times more extended 
in the black hole case. The inclination of the system has a weak 
impact on this result. The real gamma-ray photosphere would 
likely be contained between the anisotropic and the isotropic 
case i.e. 10^-10'" cm, regardless the nature of the compact object 
in Cygnus X-3. 

3. Absorption and emission in the corona 

The goal of this part is to investigate the role of a corona in 
emission and absorption of high-energy gamma rays. First, we 
estimate the contribution of non-thermal hard X-rays from the 
corona to the gamma-ray opacity ( 93. lb . Second, we model the 
emission of the corona and look whether GeV gamma rays can 
be produced and escape the system in Cygnus X-3 03.2l i. 

3.1. Absorption with non-tliermal liard X-rays 

We modeled the corona as a spherical region of radius cen- 
tered at the location of the compact object where energetic 
electrons are injected and radiate via synchrotron and inverse 
Compton scattering. For an optically thin and homogeneous 
spherical s ource, the emerging photon density h as a radial de- 
pendence riouldlll979l: iBand & Grindiavlll985h . We approxi- 
mate the emission inside the corona (p < Rc, where p is the 
distance to the center) as homogeneous and isotropic. Outside 
(p > Re), photons are assumed to escape radially, hence the flux 
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Fig. 6. Optical depth at infinity for 1 GeV photons interacting 
with non-thermal hard X-rays (10-100 keV). Gamma rays are 
injected on the axis of the accretion disk at an altitude zir = 0) 
and with an inclination t// = 30°. The corona X-ray density is 
homogeneous and isotropic in a spherical region of radius R^ = 
6 Tg (solid line), 60 rg (dashed line) and 600 (dotted-dashed 
line) around the compact object, here a 20 M© black hole (6 rg = 
1.8 X 10^ cm). Outside, the X-ray flux is radial. The total power 
injected in hard X-rays is Lc = 5 x 10^^ erg s"'. The gamma-ray 
optical depth from the disk only is shown with the red dashed 
line for comparison. 

decreases as 1/p^ and the corona appears as a point- like emitter. 
X-ray photons are injected with a power-law energy distribution 
of index -2 ranging from 10 keV to 100 keV, corresponding to a 
total power injected in the corona of Lc = 5 x 10^^ erg s ' . This 
simple spectral parametrization is a rather good representation 
for the hard X-ray tail observed in the soft state above 10 keV 
(and radio flares, see e.g jSzostek et al.ll2008i Fig. 8). AppendixlB] 
gives more details on the model adopted. 

Fig.|6]shows the contribution to the gamma-ray optical depth 
at 1 GeV integrated along the line of sight as a function of the 
location of the gamma-ray source above the accretion disk (for 
r - 0). Inside the corona, Tyy is almost constant and gamma 
rays are only marginally absorbed even if the corona is com- 
pact {Rc ~ Rin). Outside, the gamma-ray optical depth drops 
because of the decrease of the coronal photon density with dis- 
tance oc 1/p^. Gamma-ray absorption by hard X-rays dominates 
over absorption by the accretion disk photons beyond a certain 
altitude z if the corona is extended, but the optical depth is very 
small in this case (r^y «; 1) unless the luminosity in the corona 
is very high > 10^^ erg s '. However, such high flux is not 
observed in hard X-rays even during major radio flares where 
can be as high as < 5 x lO^*^ erg s"' (between 10 and 100 keV). 
This simple calculation shows that the non-thermal X-ray com- 
ponent in the soft state is not intense enough to produce an extra 
absorption feature in gamma rays. The results obtained in the 
previous section are unchanged by the presence of the corona. 

3.2. Gamma-ray emission from the corona 

In this section, we aim to model more precisely the emission 
from the corona in Cygnus X-3, and see whether it could ex- 
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Fig. 7. Simulated high-energy spectrum from the corona in 
Cygnus X-3 with BELM. The muhicolor black body spectrum 
from the accretion disk (black dotted line) is comptonized by 
a population of thermal and non-thermal electrons heated and 
accelerated in the corona up to a maximum Lorentz factor 
Tmax = 2 X 10^. The black line shows the best-fit solution 
(not corrected from absorption of soft X-rays by the ISM and 
the stellar wind) to the X-ray data (archive RXTE data, pink 
crosses from Szosteketal. 2008) in the ultra-soft state given 
in IHjalm arsdotter et al.l (l2009ft . This solution corresponds to a 
corona of radius =8x10^ cm and a total unabsorbed bolo- 
metric luminosity Lboi = 3.6 x 10""^ (at 9 kpc). This figure shows 
also the spectra for - 8x10** (green), 8x10'' (blue) and 
8 X 10"^ cm (red), assuming the same bolometric luminosity. 
The dashed lines are the intrinsic spectra (without pair produc- 
tion) and the solid lines are the absorbed spectra reaching t he 
observer. The Fermi power-law taken from lAbdo et al ] (l2009h IS 
overplotted for comparison. 

plain t he observed gamma-ray emission. In lHjalmarsdotter et alj 
(l2009l) the X-ray spectra of the ultra-soft state of Cygnus X-3 
were fitt ed with the th ermal/nonthermal Comptonization model 
EQPAIR (ICoppilll999l) . In this model the X-ray source is mod- 
eled as a spherical Comptonizing cloud in which soft photons are 
uniformly and isotropically injected with a multicolor accretion 
disk spectrum. The electron distribution of the Comptonizing 
electrons is hybrid: at low energy the distribution of electrons 
is Maxwelllian with a non-thermal quasi-powerlaw extension at 
high energies. The Compton cloud is powered through two chan- 
nels: 1) direct heating of the thermal electrons, 2) injection of 
high-energy electrons with a power-law distribution. The steady 
state particle and photon distributions in the cloud are then com- 
puted according to energy and electron/positron pair balance, as 
a function of the power input in the corona, as well as other pa- 
rameters, such as the temperature of the soft seed photons, size 
of the emitting region or Thomson optical depth t of the plasma. 
The calculation is performed in the one-zone approximation and 
radiative transfer is dealt with using an escape probability for- 
malism. 

From this analy si s of t he X-ray spectra of Cygnus X-3, 
iHialmarsdotter et all (l2009b infer that in the ultra-soft state, 
about half of the coronal power is provided to the electrons in 
the form of non-thermal electrons injected with a slope Finj - 
2.1. In principle, these non-thermal electrons could also be re- 
sponsible for the gamma-ray emission detected by Fermi. This 



would require that the coronal electrons are accelerated up to 
a Tmax > 10'*. In the EQPAIR model the grids used to compute 
the electron and photon energy distributions are fixed and can- 
not be exten ded to an energy high enough to simulate the Fermi 
emission. In IHialmarsdotter et al.l (12009.) the maximum Lorentz 
factor of the injected electrons is ymax = 10^, clearly too low 
to produce the gamma-ray emission. However the exact value 
of Tmax is not constrained by the X-ray data and could be much 
larger 

In the following, we use the BELM code (iBelmont et alJ2008h 
to estimate the gamma-ray flux predicted by this model. The 
BELM code is very similar to EQPAIR (i.e. based on the resolution 
of the kinetic equations in the one-zone, isotropic approxima- 
tion) but offers more flexibility regarding the maximum energy 
of the electrons. BELM can also take into account synchrotron 
emission and self-absorption (in addition to Compton, electron- 
proton Bremsstrahlung emission, pair production/annihilation 
and the internal electromagnetic c ascade of pairs). We use ex- 
actly the same parameters found by IHialmarsdotter et al.l(l2009h 
in their best fit model of the ultra-soft state except that the max- 
imum Lorentz factor is now fixed at y^ax = 2 x 10^ in order 
to produce gamma-ray emission up to 100 GeV. The total power 
injected as non-thermal electrons is set to Lnth = 3 x 10^^ erg s"' . 

We find that high-energy gamma rays are mostly produced 
by inverse Compton scattering. However, gamma rays suf- 
fer from a strong absorption via pair production with non- 
comptonized disk photons (see Fig. |7]l. The gamma-ray flux 
cannot be explained. The radiation from the cascade of pairs 
in the corona is not sufficient to increase the gamma-ray flux 
to an observable level. Pair production can be reduced if the 
corona is more extended. The best fit solution to the X-ray 
data in the ultra-soft state corresponds to a corona of radius 
/Jc = 8 X 10^ cm ^ 30 rg (for a 20 Mq black hole). By keeping the 
total bolometric luminosity constant (Lboi = 3.6 x lO-''^ erg s"'), 
the level of emission and the spectrum in the Fermi energy band 
can be reproduced if the radius of the corona R^ > 10'*' cm 
(Fig.EJ. 

These results are consistent with the detailed calculation of 
gamma-ray absorption above the disk (i.e. inhomogeneous in the 
full isotropization limit, see Sect. 12.4b . Some spectral changes 
appear in X-rays because of the reduction in the compactness 
of the corona. Fig. |7] shows it will be diflicult to reconcile the 
flat X-ray spectrum at 100 keV seen in very soft states (pro- 
duced by electrons when the compactness is high) with signif- 
icant gamma-ray emission (which requires low compactness to 
avoid gamma-ray attenuation by pair production). Note also that 
the spectrum could be changed due to anisotropic effects and in- 
homogeneities not taken into account in this model. It is difficult 
to imagine a very extended corona since most of the gravitational 
potential energy (90%) remains in the region within 30 r^ around 
the compact object. 

The simulations in Fig. |7] assumes no magnetic field. For 
magnetic energy density at equipartition with the radiation in 
the corona, the magnetic field is about Z?eq ~ 5 x lO*" G (for 
/?c = 8 X 10^ cm). Synchrotron radiation becomes the dominant 
processes at high energies, but the shape of the escaping high- 
energy spectrum changes only slightly. The energy of the absorp- 
tion cut-off remains magnetic field independent, hence the con- 
straints on the size of the corona are unchanged. In addition, the 
high-energy cut-off" of the intrinsic synchrotron spectrum can- 
not exceed WRC' /gp ^ 70 MeV (wh ere is the fine structure 
constant, see e.g. lGuilbert et al.lll983h due to synchrotron losses, 
hence too low to account for the high-energy emission. This limit 
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is calculated assuming no relativistic Doppler boosting and that 
the cooling timescale via synchrotron emission tsyn equals the ac- 
celeration timescale of particles in the corona which cannot be 
shorter than the Larmor timescale ti^, assumption usually made 
e.g. in diffuse shock acceleration. This balance between acceler- 
ation and cooling implies also tha t it is very difficult to acceler- 
ate particles at high energies (e.g. 'Hillas 1984: 'Ah aronian et alJ 
I2OO2.) . For the equipartition magnetic field strength, the condi- 
tion Uyn > ?acc gives y^.^^ < (9m^c'*/4e^Beq)'^^ = 5 x 10"*, allow- 
ing gamma-ray emission to a few tens of GeV only. To accelerate 
particles to y^ax =2x10^, the magnetic field should not exceed 
~ 3 X 10^ G. These estimates are based on the optimistic assump- 
tion that the acceleration efficiency defined as // = E^jeBc equals 
1 (where is the synchrotron energy losses and e is the charge 
of the electron). The constraint on the magnetic field strength in 
the corona could be even more stringent than what is given here. 
The magnetic field strength at equipartition is also sufficiently 
high to quench the development of pair cascade. 

We conclude that the gamma-ray emission measured by 
Fermi has probably not a coronal origin. 



likely to be produced by energetic electron-positron pairs pos- 
sibly located in the jet, upscattering stellar radiation via inverse 
Compton scattering. In this case, the Fermi spectrum is well fit- 
ted if pairs are injected with a power-law energy distribution of 
index 4.4. The power-law cannot extend below y^in > 50 (as- 
suming no Doppler boosting) or the observed hard X-ray flux 
below 100 keV would be overestimated. Alternatively, the in- 
verse Compton emission from the jet could also be responsible 
for the hard X-ray emission (instead of the corona). This would 
require a spectral break in the electron energy distribution be- 
low ybrk ^10^ with a harder power-law of index close to 3 with 
Tmin S 10 (in order to obtain a flat spectrum above 10 keV). 
The emerging picture is that high-energy gamma rays should be 
emi tted far from the a ccretion disk, probably in the relativistic 
jet dPubus etalj|2010l) . 
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4. Conclusion 

High-energy gamma-ray emission in Cygnus X-3 seems to be 
related to the soft X-ray state, i.e. when the X-ray spectrum is 
dominated by a bright accretion disk component. Gamma rays 
should be significantly absorbed by soft X-rays emitted by the 
inner (hottest) parts of the accretion disk, unless the gamma-ray 
source lies far enough from the disk. The apparent lack of ab- 
sorption feature in the spectrum measured by Fermi indicates 
that high-energy gamma rays should be located at least ~ 10**- 
10^° cm away from the accreting star, depending on the inclina- 
tion of the system. These conclusions are not affected by absorp- 
tion of X-rays in the stellar wind, since this process becomes sig- 
nificant at larger distances > 10'" cm. Thomson scattering with 
free electrons in the wind could redistribute and isotropize the 
X-ray emission above the disk and affect the gamma-ray optical 
depth map above the accretion disk. In the limit where X-rays 
from the disk are fully isotropized, the gamma-ray photosphere 
is quasi-spherical and is located at about 10"^ cm from the com- 
pact object. A more detailed treatment of this effect would re- 
quire for instance Monte Carlo techniques to follow the multiple 
elastic scattering of X-rays in the wind. 

In addition to the bright black body emission from the ac- 
cretion disk, the X-ray spectrum contains a non-thermal tail in 
hard X-rays probably related to a corona above the disk. We 
modeled the corona as a spherical, isotropic and homogeneous 
cloud around the compact object. We showed that gamma-ray 
absorption with non-thermal hard X-rays from the corona is 
too small compared with the contribution from the disk. In this 
article, we present also a more precise modeling of the non- 
thermal radiation from the corona in Cygnus X-3 with BELM. 
This code enables to extend the energy of the electrons in 
the hig h-energy domain. U s ing th e best fit solution to X-ray 
data of iHialmarsdotter et al.l (l2009h . we found that gamma rays 
produced in the corona suffer from strong absorption by non- 
comptonized photons from the disk before they escape to the ob- 
server The gamma-ray flux seen by Fermi could be reproduced 
if the corona is unrealistically extended (/?c S 10'° cm). This 
study does not favor a coronal origin of the observed gamma-ray 
emission. 

At distances comparable to the orbital separation, the Wolf- 
Rayet star photon density dominates over the accretion disk X- 
ray density. The gamma-ray emission and the modulation are 



Appendix A: Calculation of the gamma-ray optical 
depth in the radiation field of a standard 
accretion disk 

We aim to give here the detailed calculation of the optical depth 
Tyy for a gamma-ray photon of energy ei propagating above an 
accretion disk. The disk is optically thick, geometrically thin and 
flat with an inner radius and external radius /?out as depicted 
in Fig. lA.ll The compact object lies at the center of the disk 
at the point O. We consider a point-like source of high-energy 
gamma rays located at a radial distance r, at an altitude z above 
the accretion disk and with an angle 0, with the x-axis in the 
plane of the disk. A distant observer sees the disk with an angle 
i/r, this direction is set in the (xOz) plane. 



A. 1 . Total gamma-ray optical depth 

The total gamma-ray optical depth Tyy integrated over the high- 
energy photon path length / and the target photon density from 
the disk dn/dedD. (in ph cm"^ erg"' sr~') yields 

r""-" r+'" dn 

Tyy- I I I (I - COS 60) CTyyded^dl, (A.l) 

Jo Jo Je,„i„ dedD. 

where e is the disk photon energy, 0o the angle between the 
two photons, o-T,y the p air production cross section (see e.g. 
Gou ld & Schrededl 19671) and Qdisk the total solid angle covered 
by the accretion disk seen from the gamma ray. Eq. (IA.lt must 
also satisfy the kinematic condition 

2my 

e > emin = — -; --t (A.2) 

ei (1 - cos 60) 

where m^ is the rest mass of the electron (or positron), otherwise 
no pairs are produced i.e. Tyy - 0. 



A.2. Geometrical aspects 

In order to derive explicitly the gamma-ray optical depth in 
Eq. (lA.lb . a few geometrical relations should be established. 
First, we define the following unit vectors in cartesian coordi- 
nates as (see Fig. lA.ll) 



Cobs = (sin ij/, 0, cos i//) 



(A.3) 
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Perspective view 



Top view 



Observer 





Observer 



Fig. A.l. Diagram of a flat, geometricany fliin accretion disk formed around the compact object in the center (point O). A source of 
energetic gamma rays Hes above the disk seen at an angle by a distant observer. Gamma rays and photons from the unit surface of 
the disk dS (centered around point M) interact at the point P with an angle Oq. Point Q is the projection of point P in the disk plane. 
All the other geometrical quantities referred in the text for the calculation of the gamma-ray optical depth are also shown here. 



along the direction joining the source and the observer, 

Cr - (cos (p, sin 4>, 0) (A.4) 

along the OM direction, 

Cp = (sin 0COS (pp, sin 6 sin (pp, cos 9j (A. 5) 
along the OP direction and 

Cdisk = (sin a cos oj, sin a sin oj, cos a) (A. 6) 
in the MP direction. Elementary trigonometry gives 

cos0^- -, sin6»=-!^ —, (A.7) 



P 



P 



where x - r cos <ps and y = r sin 0s are the coordinates of the 
source in the plane of the disk. Combining the expressions given 
in Eq. (IA.7b . we obtain 



- + + + 21 (x sin i/r -H z cos i/r) . (A.8) 

The distance D = MP between the element surface of the disk 
dS and the point where both photons collide is 



= + p-2Rpcosx, 
with 

cosx = Cp • Cr = sin COS ipp cos <p + sin 6? sin <pp sin <p, (A. 10) 
where 



X + Isini// y 
cos <j>n : , sm <f>n - 



psinO 



psinO 



(A. 11) 



The angle a between the direction MP and the axis of the disk z 
introduced in Eq. (IA.6b is defined such as 

pcose 5 (D^-p^cos^of^ 

cosQ- , sina = — = . (A. 12) 

D D D 



The cosine of the scattering angle 0o between the high-energy 
gamma ray and the photons from the element surface of the disk 
is 

cos 00 = Cobs • Cdisk - sin if/ sin a cos u + cos ijjcosa, (A. 1 3) 
where 

X + ls\n\j/ - R cos (p 



cos 0) 



{D^ - p2 cos 



,2 nr,c2flll/2 



(A. 14) 



The element solid angle dD. covered by the surface dS - RdRd(p, 
as seen by a gamma-ray photon along its way to the observer is 
(see Fig. lA.lb 

dD. = ^2 = ^-—^RdRdcf>, (A. 15) 

Cz is the unit vector along the z-axis. 

A.3. Soft photon density of the disk 

In the standard model, the accretion disk is formed by concentric 
annuU i n thermal equ ilibrium emitting a black body spectrum 
(see e.g. lPringlel[T98ll) . Hence, a surface element of the disk dS 
generates the density 



dn 



(A.9) dedQ. h^c^ exp [e/kT(R)] - 1 



(A. 16) 



with h, the Planck constant, k the Boltzmann constant and T{R) 
is the eff'ective temperature of the disk at the distance R from 
the center. Using the temperature profile given in Eq. ([1), the 
spectrum integrated over the disk is 



dn _ r""-" 
de Jo 



dn 
dedQ. 



dQ, 



(A. 17) 



where dO. is given in Eq. ( lA.lSb . In the case of an observer lo- 
cated at a distance L » so that D ^ p = L and k if/, dQ. 
simplifies into dO. ~ cos t//RdRd(p/L^. The isotropic apparent 
flux of the disk seen by the observer eL^ (in erg s ') is then 



eLf = AncL^e^'^ = Stt^c cos i/r 
de 



dn 
dedQ. 



RdR. 



(A. 18) 
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In the limit Ri^ <K 7? ^ Rout, the flux can be expressed as 



37r2/!c2/3 



■ cos i/r M^/^M^fe^/a^ 



a power-law of index 4/3 in energy. G is the gravitational con- 
stant. ^ and r are respectively the Riemann and the gamma func- 
tions, withf(8/3)« 1.2842 and r (8/3) ~ 1.5046. 



Appendix B: Calculation of the gamma-ray optical 
depth with non-thermal hard X-rays from the 
corona 

B.1. Photon density of the corona 

Hard X-rays are assumed to originate from the corona, a spheri- 
cal region of radius R^ (see Fig. IB. II ). The density of non-thermal 
X-rays follows 



dn 
ae 



(B.l) 



where K^^ is a normalization constant and p the spectral index. 
For simplicity, we assume that X-rays are emitted radially out- 
side the corona (if p > R^). If is the total luminosity of the 
corona, we have 



r* dn 
\ e — de. 
de 



(B.l) 



Hence the density of non-thermal X-rays outside the corona is 

\2 



driest . . / , _„ 

where (if p i^2) 

(2-p)U 



47Tc{e-/^'-e:"^')Rl' 



(B.3) 



(B.4) 



If X-ray photons are homogeneous and isotropic inside the 
corona (p < Rc), then 



dnini drif. 



de de 



(p^Rc)^ Koe-P. 



B.2. Optical depth 



Gamma rays are injected on the axis of the accretion disk at an 
altitude z. If z < R^, gamma rays first cross the corona over a 
distance lo given by 



/() = \^rI-£-{\ -cosV)]'^^ -zcos(/r 



(B.6) 



and escape at the point N (see Fig. IB. 11 1. We average over an 
isotropic distribution of pitch angle between both photons 9q. 
The integral over the length path in the full expression of the total 
optical depth (Eq. lA.lb is trivial since there is an invariance in I 
for homogeneous target photon density. The gamma-ray opacity 
inside the corona is 



Tint (ei 



de 



(1 - cos 9q) (Tyy sin OodOode. (B.7) 



Observer 




Fig. B.l. The corona is approximated as an optically thin spher- 
ical region of radius R^ (violet region) centered at the com- 
pact object (point O), and emit hard X-rays. Inside the corona 
(p < Rc), the density of hard X-rays is uniform and isotropic. 
Outside (p > Rc), photons propagate outward radially. Gamma 
rays are injected on the disk axis with an inclination ifr towards 
the observer 



Outside the corona. X-ray photons are anisotropic and not ho- 
mogeneous. The angle of interaction is related to the length path 
/ via 

cos 6*0 = Cobs ■ Cp 

/(/o + OsiniA „ z+(;o + OcosiA\ 
= (sm (/r, 0, cos (/r) ■ ,0, 



Iq + I + z cos l// 

p 



(B.8) 



where Iq is given in Eq. (IB. 61) . I is the distance NP, and p is the 
distance to the interaction point from the center OP such that 



p2 = Ho + l)^ + +2z (lo + I) cos i/r. 
The gamma-ray optical depth is then 

dricn 



(B.9) 



T ext (ei 



Jl„ Je_ 



de 



(p) (1 - cos 9o) o-yydedl. (B.IO) 



The total optical depth is Tyy - Tint + Text 
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